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ABSTRACT: We study an information-theoretic measure of uncertainty for quantum 
systems. It is the Shannon information / of the phase space probability distribution {z\p\z), 
where \z) are coherent states, and p is the density matrix. As shown by Lieb and Wehrl, 
I > 1, and this bound represents a strengthened version of the uncertainty principle. 
^ ■ For a harmonic oscillator in a thermal state, I coincides with von Neumann entropy, 
— Tr(plnp), in the high-temperature regime, but unlike entropy, it is non-zero (and equal 
to the Lieb- Wehrl bound) at zero temperature. It therefore supplies a non-trivial measure 
of uncertainty due to both quantum and thermal fluctuations. We study I as a function of 
time for a class of non-equilibrium quantum systems consisting of a distinguished system 
coupled to a heat bath. We derive an evolution equation for I. For the harmonic oscillator, 
in the Fokker-Planck regime, we show that / increases monotonically, if the width of the 
coherent states is chosen to be the same as the width of the harmonic oscillator ground 
state. For other choices of the width, and for more general Hamiltonians, I settles down 
to monotonic increase in the long run, but may suffer an initial decrease for certain initial 
states that undergo "reassembly" (the opposite of quantum spreading) . Our main result is 
to prove, for linear systems, that / at each moment of time has a lower bound J™ 7 ^ over 
all possible initial states. This bound is a generalization of the uncertainty principle to 
include thermal fluctuations in non-equilibrium systems, and represents the least amount 
of uncertainty the system must suffer after evolution in the presence of an environment for 
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time t. 7™ n is an envelope, equal, for each time t, to the time evolution of I for a certain 
initial state, which we calculate to be a squeezed coherent state. Jj nin coincides with 
the Lieb-Wehrl bound in the absence of an environment, and is related to von Neumann 
entropy in the long-time limit. The form of I™ m indicates that the thermal fluctations 
become comparable with the quantum fluctuations on a timescale equal to the decoherence 
timescale, in agreement with earlier work of Hu and Zhang. Our results are also related 
to those of Zurek, Habib and Paz, who looked for the set of initial states generating the 
least amount of von Neumann entropy after a fixed period of non-unitary evolution. 
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I. INTRODUCTION 



One of the most important features of quantum mechanics is the uncertainty principle, 

AxAp > ^ (1.1) 

Although frequently interpreted as a statement about the precision of measurements, it 
may also be taken to mean that there is intrinsic uncertainty in any phase space description 
of quantum systems. This uncertainty may be especially significant for systems in certain 
states, such as the ground state. However, in many quantum systems of interest there is 
additional uncertainty due to thermal fluctuations, and moreover, there may be regimes 
in which the thermal fluctuations dominate. A number of questions then naturally arise: 
Is there a useful measure of uncertainty due to both quantum and thermal fluctuations? 
And, if so, what is the lower bound on this uncertainty, analogous to (1.1)? What are the 
regimes in which each type of fluctuations dominate? This paper addresses these questions. 

Apart from being of interest in their own right, there are a number of specific moti- 
vations for studying these issues. The principal one concerns the general question of the 
emergence of classical behaviour in quantum systems. Understanding this issue is one of 
the main aims of the decoherent histories approach to quantum mechanics [1,2,3,4]. There 
(and in other approaches [5,6,7,8,9]), the process of decoherence is held to play an essential 
role. This process typically occurs as a result of interaction of the system under scrutiny 
with a wider environment. But this same interaction also leads to essentially random dis- 
turbances of the system, driving it off its classical path. The probabilities for histories 
are typically found to be peaked about classical histories, with some width determined by 
quantum effects and broadened by thermal fluctuations induced by interaction with the 
environment [1]. It therefore becomes important to gain a quantitative understanding of 
both types of fluctuations, and to find the regimes in which each are important. 
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In this paper we will explore an information-theoretic measure of uncertainty due to 
both quantum and thermal effects, suitable for the non-equilibrium quantum systems used 
in decoherence models. 

We begin in Section II by describing the necessary background. We first review some 
aspects of information theory. We then introduce a quantum-mechanical phase space 
distribution. It is the distribution 

fi(p,q) = (z\p\z) (1.2) 

where p is the density matrix of the system, and \z) are the coherent states. Our chosen 
measure of uncertainty is the Shannon information / of this distribution, 

I = ~ J 7^ /*(P> q) In/xO, q) (1.3) 

As we shall explain, the uncertainty principle manifests itself through the inequality, 

I > 1 (1.4) 

with equality if and only if p is a coherent state [10,11]. Our main aim is to generalize 
(1.4) to include the effects of thermal fluctuations in non-equilibrium systems. 

In Section III we study the properties of I for a simple equilibrium system - the 
harmonic oscillator in a thermal state. This simple example clearly illustrates how / 
supplies a useful measure of both thermal and quantum fluctuations. We then go on, in 
Section IV, to consider non-equilibrium systems, the main topic of this paper. We describe 
an important class of non-equilibrium systems consisting of a distinguished system coupled 
to a heat bath (often referred to as open quantum systems). 

In Section V, we discuss the time-evolution of I for non-equilibrium systems. We show 
that It generally settles down to monotonic increase. There is, however, the possibility 
of an initial period of decrease for specially chosen initial states which reassemble (the 
opposite of wavepacket spreading). 
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In Section VI, we describe our main result. This is the demonstration that It has a non- 
trivial lower bound, the generalization of the Lieb-Wehrl result (1.4) to include thermal 
fluctuations in non-equilibrium systems. The function I™ m bounding It from below is 
generally not the time evolution of / for some particular initial state, but is an envelope. 
The initial state which achieves I™- 171 at time t (but generally not at any other time) is 
a squeezed coherent state, with a specific value for the squeezing factor depending on t. 
jrnin j g a measure Q f the least amount of quantum and thermal noise the system must 
suffer after non-unitary evolution for time t. The bound reduces to the Lieb-Wehrl bound 
in the absence of an environment. 

As we shall explain, there are three contributions to the uncertainty: 

(1) There is the uncertainty intrinsic to quantum mechanics, expressed through the uncer- 
tainty principle, (1.1). This is not dependent on the dynamics. It is this uncertainty 
that is referred to by the expression "quantum fluctuations" . 

(2) There is uncertainty that arises due to the spreading or reassembly (the reverse of 
spreading) of the wave packet. This effect depends on the dynamics, and because 
quantum mechanics is time-symmetric, it may increase or decrease the uncertainty. 

(3) There is the uncertainty due to the coupling to a thermal environment. This has two 
components: dissipation and diffusion (the latter being responsible for the process of 
decoherence) . This generally tends to increase the uncertainty as time evolves. 

The point is that the lower bound, IJ nm : includes the effects (1) and (3), but avoids (2). 

Finally, in Section VII, we summarize and discuss our results. We compare our results 
with calculations of Hu and Zhang [12], who calculated the time evolution of the usual 
uncertainty function for a particular initial state, and determined the timescale on which 
the thermal fluctuations catch up with the quantum fluctuations. We also compare with 
the results of Zurek, Habib and Paz [7,8], who looked for the set of initial states which 
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generate the smallest amount of von Neumann entropy after a fixed period of non-unitary 
evolution. 



II. BACKGROUND 

We now review the necessary background. 

11(A). Information Theory 

Suppose one has a set of probabilities pi for a data set S consisting of discrete set of 
alternatives labeled by i, i = 1,2- • -N. One has < pi < 1 and J2iPi = 1- Then the 
Shannon information of the data set is defined to be 

N 

I(S) = -5>ln Pi (2.1) 

i=l 

Here, In is the logarithm to base e. I(S) satisfies the inequalities 

0</(S)<lniV (2.2) 

It reaches its minimum if and only if p^ = 1, for one particular value of z, and so pi = 
for all the other values. It reaches its maximum when Pi = for all i. The information 
of a probability distribution is therefore a measure of how strongly peaked it is about 
a given alternative. For this reason, I(S) is sometimes referred to as uncertainty, being 
large for spread out distributions and small for concentrated ones. This nomenclature is 
appropriate for purposes of this paper. The expression (2.1) is also often referred to as the 
entropy of the distribution, but we will not do so here, reserving the word entropy for the 
von Neumann entropy of quantum statistical mechanics (discussed in later sections). 
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In a similar manner for continuous distributions, let X be a random variable with 
probability density p(x). Then J dx p(x) = 1. The information of X is defined to be 

I{X) = - J dx p(x) lnp(x) (2.3) 

Unlike the discrete case, I(X) is no longer positive, since p(x) is not a probability, but a 
probability density, so may be greater than 1. However, it retains its utility as a measure 
of uncertainty. This is exemplified by a Gaussian distribution of variance Ax, 

P(*) = " r ex P f-fc^T) (2-4) 

(2tt(Ax)2)2 V 2(Ax)2 ; 

It has information 

i 

I{X) = In (27re(Aa0 2 ) 2 (2.5) 

From this we see that /(X) is unbounded from below, and indeed, approaches — oo as 

Ax — > and p(x) approaches a delta-function. I(X) is also unbounded from above, as 

may be seen by taking the width Ax to be very large. However, if the variance is fixed, then 

a straightforward variational calculation shows that I(X) is maximized by the Gaussian 

distribution (2.4). We therefore have the important inequality, 

i 

I{X) < In (^27re(Ax) 2 ) 2 (2.6) 

The generalization to probability distributions of more than one variable is straight- 
forward. For example, one has, 

I(X, Y) = - J dxdy p(x, y) \np{x, y) (2.7) 

and it is easy to show that 

I(X,Y) < I(X) + I(Y) (2.8) 

where I(X) is the information of the distribution J dyp(x, y), and similarly for I(Y). We 
also record another useful result. Let f(x),g(x) > and let J dxg(x) = 1. Then 

J dxf{x)g{x)^j\n(^j dyf(y)g(y^j > - J dxf(x)g(x)lnf(x) (2.9) 



This is essentially due to the convexity of the function xlnx, and also holds in the discrete 
case. Further details on information theory may be found in the literature [13]. 



11(B). Phase Space Distributions in Quantum Mechanics 

As stated above, our work is partly aimed at discussing the emergence of classical 
behaviour. In this connection, it is often useful to introduce quantum-mechanical phase 
space distributions. There are a variety of phase space distributions that may be employed 
in quantum mechanics [14]. In this paper we shall focus on the function, 

fi(p,q) = (z\p\z) (2.10) 

where 

(x\z) = (x\p, q) = exp [~^- + 

are the coherent states, with a p a q = ^h. We find it useful to work with units with 
dimension, and for this reason it is necessary to introduce the parameter a q into the 
coherent state wave functions. The function fi(p, q) is normalized according to 

f dpdq 

It is readily shown that p(p, q) is also equal to 

p(p, q) = 2 J dp'dq' exp " ^T^j W P& > ^ C 2 - 11 ) 

where Wp(p, q) is the Wigner function of p, defined by [14], 

W p (p, v) = ^-JdZ e-irt p{q + ±Z,q- (2.12) 

The distribution p(p, q) is therefore a Wigner function, smeared over an K sized region of 
phase space. This smearing renders the distribution function positive, even though the 
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Wigner function is not in general [15]. The distribution (2.11) is sometimes known as the 
Husimi distribution [16], and has appeared frequently in discussions of the Wigner function 
(e.g, Refs.[l, 15,17]). 

The utility of the distribution function ft(p, q) will become apparent as we expose some 
of its properties. We remark, however, that \x is of the form, 



where P z = \z)(z\ is a coherent state projector (actually only an approximate projector due 
to the over-completeness of the coherent states), fi(p, q) therefore has the interpretation as 
the probability of a simultaneous but approximate sampling of position and momentum. 
Moreover, it may be shown that that by taking suitably weighted sums over p and q of 
(2.13), an object of the form 



may be obtained, where Px(t) denotes an imprecise position sampling at time t. Eq.(2.14) 
is the probability for the history characterized by the initial state p, and samplings of 
position at times t\ and The distribution p(j>, q) is therefore closely connected with the 
the decoherent histories approach to quantum mechanics, which focuses on objects of the 
form (2.14). In particular, it may be shown that the degree to which expressions of the 
form (2.14) are peaked about classical paths is limited by the degree of peaking of fi{p, q) 
in phase space. This is discussed in another paper [18]. 



fi{p, q) = Tr[P z p] 



(2.13) 



P (x2,t2,x 1 ,t 1 ) = Tr[Px 2 {t 2 )Px 1 (tl)pPx 1 {ti)} 



(2.14) 



11(C). An Information-Theoretic Measure of Uncertainty 



We are interested in the extent to which fi(p, q) is peaked about some region of phase 
space. As we have discussed, the Shannon information is a natural measure of the extent 
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to which a probability distribution is peaked. We shall therefore take as our measure of 
uncertainty, the information 

I(P, Q) = -Jt^ ^ q ) i^fo q ) (2.15) 

The uncertainty principle strongly suggests that a genuine phase probability distribu- 
tion in quantum mechanics cannot be arbitrarily peaked about a point in phase space. 
We therefore expect the information (2.15) to possess a lower bound. Furthermore, since 
coherent states are normally regarded as the states most concentrated in phase space, we 
expect the lower bound to be the value of / on a coherent state. It turns out that both 
of these expectations are true. It was conjectured by Wehrl [10], and proved by Lieb [11], 
that 

I(P,Q) > 1 (2.16) 

with equality if and only if p is the density matrix of a coherent state, \z )(z \ . 

The inequality (2.16) may be related to the usual uncertainty principle, (1.1). One has 
the inequalities, 



ln(|A^A^) > I(Q) + I(P) 



> I(P,Q) (2.17) 

The second inequality is an elementary property of information, (2.8); the first is the 
inequality (2.6) applied to each of the marginal distributions for p and q, where A^q and 
A^p are the variances of the distribution /j,(p, q) (the difference by a factor of 27rh is due 
to our choice of phase space measure). These variances are, however, not the quantum- 
mechanical variances, since they include the variances of the coherent states. Indeed, one 
has 

(A^) 2 = (A^) 2 + a g 2 (2.18) 

{A^p) 2 = {A p pf + al (2.19) 
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where A p denotes the quantum-mechanical variance. Now (2.16)-(2.19) together imply 
that 

[(A p qf + a 2 ) ((A p p) 2 + a 2 ) > h 2 (2.20) 

Now note that the width a q in the coherent state is so far arbitary. Minimizing (2.20) 
over a q (and recalling that a q a p = \h), we thus obtain the standard uncertainty relations, 
(1.1). An alternative method of connecting the standard uncertainty relations with (2.16) 
has been given by Grabowski [19]. 

Suppose now we have a state which is genuinely mixed. It may therefore be written, 

P = ^Pn\n){n\ (2.21) 

n 

for some basis of states |n), and where p n < 1. One has 

n 

The information of (2.22) will always satisfy (2.16), but this will be a very low lower bound 
for a mixed state. However, from the inequality (2.9), one has 

n 

= Pn Inpn 

n 

= -Tr(plnp) = S[p] (2.23) 

That is, / is bounded from below by the von Neumann entropy, S[p\. As we shall see in 
the following section, this inequality can be close to equality in the regime where thermal 
fluctuations are large. This close connection with von Neumann entropy is one of the 
virtues of our chosen measure of uncertainty, over other measures one might contemplate 
(e.g., the usual uncertainty function, U = (A p q) 2 (App) 2 ). 

From the above, we therefore see that / is a useful of measure of both quantum and 
thermal fluctuations. It possesses a lower bound expressing the effect of quantum fluctu- 
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ations, and is closely connected to entropy, which in turn is a measure of thermal fluctu- 
ations. In the following sections we will explore the further properties of /, especially for 
non-equilibrium systems. 



III. FLUCTUATIONS AT THERMAL EQUILIBRIUM 

To see some of the features of I more clearly, consider the equilibrium case. Let the 
density matrix be thermal, p = e~^ H , where Z = Tr(e~P^) is the partition function, 
and /3 = 1/kT. One has 

n 

where \n) are a set of energy eigenstates with eigenvalues E n . For simplicity, we restrict 
attention to the simple harmonic oscillator, for which, 

and so E n = ?ko(n + 3), and 

I | In 

|<z|n)| 2 =£L-e-W a (3.3) 

lb . 

Here, z = ^ {q/&q + ip/ffp) > where a q a p = ^h, and we have made the choice a q = 
(h/2Muj)2 . See Ref.[20] for details about the coherent states. One thus has 

H(q,p) = (z\p\z) = (1 - e~^) exp (-(1 - e-^ huJ )\z\ 2 ^j (3.4) 

The information (2.15) may then be computed explicitly. It is, 

I = 1 - In (l - e"^) (3.5) 

Eq.(3.5) is exactly the sort of result one would expect. As the temperature goes to zero, 
(3 — > oo, and the uncertainty reduces to the Lieb-Wehrl result, (2.16), expressing purely 
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quantum fluctuations. But for non-zero temperature, the uncertainty is larger, tending 
to the value —ha((3Ku), as the temperature goes to infinity. This limit expresses purely 
thermal fluctuations. For more general Hamiltonians, we expect the information / of the 
equilibrium thermal state to behave similarly (although we have not been able to derive 
its explicit form). 

It is of interest to compare (3.5) with the entropy, 

S = -Tr(plnp) (3.6) 
The partition function is readily shown to be, 



Z = \ (3.7) 



and the entropy is 



S = -/3^(lnZ) + lnZ 

= - In ( 2 sinh(^/iw) J + \(3hu coth{]-(3hu) (3.8) 
For large temperatures (small /3), 

-\n{j3fiuj) (3.9) 

S therefore coincides with / in the high-temperature limit, On the other hand, S — > as 
the temperature goes to zero, whilst / goes to a non-trivial lower bound. 

We therefore see that I is a useful measure of uncertainty, in both the quantum and 
thermal regimes. Entropy, by contrast, supplies a measure of uncertainty due only to 
thermal fluctuations. It is therefore good in the thermal regime, but in the quantum 
regime, it underestimates the intrinsic quantum uncertainty since it goes to zero for pure 
states. 

It is also useful to compare this measure of uncertainty with the more standard measure, 

U = (A p qf (A p pf (3.10) 
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Here, (A p q) 2 is computed using (q 2 ) = Tr(g 2 p), etc. One readily finds that 

(A p p) 2 = u 2 (A p qf (3.11) 

Now Eq.(3.4) is product of Gaussians in p and q, with variances A p q, A^p, say. The 
information of such a distribution may be written, 

I = In (|A M(7 A M p) (3.12) 

As in (2.18), (2.19), the variances of q and p in (3.12) are not the same as the quantum- 
mechanical variances, because they also include the variances of the coherent state: 

(A^) 2 = (A p q) 2 + A (3.13) 
(A p p) 2 = (A p pf + ^ (3.14) 

Inserting these in (3.12) and using (3.11), one obtains, 

(3.15) 

This shows that, in this simple case, there is a complete equivalence between U and I 
as measures of uncertainty. We do not expect this equivalence to hold more generally, 
however. 

Finally, we note that an information-theoretic uncertainty relation including the effects 
of thermal fluctuations at thermal equilibrium has been derived by Abe and Suzuki [21], 
using thermofield dynamics. Their information-theoretic measure is different to the one 
used here. 



In 



IV. NON-EQUILIBRIUM SYSTEMS 



Consider now the case of non-equilibrium systems, the main topic of this paper. An 
important class of such systems in the present context are those in which the total system 
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naturally decomposes into a distinguished system, S say, and the rest, summarily referred 
to as the environment. S is then often referred to as an open quantum system. One is 
interested only in the behaviour of S, and not in the detailed behaviour of the environment. 
The distinguished system is most completely described by the reduced density matrix, p, 
obtained by tracing out over the environment. The environment leaves its mark, however, 
in that the effective evolution of the reduced density matrix alone is non-unitary. 

A useful model of the type described above consists of a particle moving in one- 
dimension in a potential V(x), linearly coupled to a bath of harmonic oscillators in a 
thermal state. The environment is characterized by a temperature T and a dissipation 
coefficient 7. This model has been the subject of many papers, so we will give only the 
briefest of accounts here (for further details, see Refs. [22,23,24,25,26,27]). 

After tracing out the environment, the reduced density matrix p of the distinguished 
system evolves non-unitarily, according to the relation 

Pt{x, V) = J dx d V0 J (x, y, t\x , y , 0) po(x , y ) (4.1) 

Here, J is the reduced density matrix propagator. It is given by the path integral expres- 
sion, 

J(x f ,y f ,t\x ,y ,0) = J VxVy exp (j.S[x] - % -S[y] + % -W[x,y^J (4.2) 



where 



^Mx 2 - V(x) 



(4.3) 



S[x] = J dt 

and W[x(t),y(t)] is the Feynman- Vernon influence functional phase, 

W[x(t),y(t)} = - f ds I" ds'[x{s)-y{s)\ r,(s - s') [x(s') + y(s')} 
JO JO 

+ i f ds f ds[x{s)-y{s)}u(s- s)[x{s)-y(s)} (4.4) 

Jo Jo 

The explicit forms of the non-local kernels 7] and v may be found in Refs. [26,23]. We have 
assumed, as is typical in these models, that the initial density matrix of the total system 
is simply a product of the initial system and environment density matrices. 
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Considerable simplifications occur in a purely ohmic environment at high temperature. 
Take a regularized ohmic environment with cutoff frequency A having the spectral density 

C{u) = 2 _M^ e -^/A\ (4.5) 

7T 

In the Fokker-Planck limit (see Ref.[26]), one first takes the high temperature limit Ti/kT <C 
A -1 and then lets the cutoff go to infinity, A — > oo. One finds 

r}{s - s) = M 1 S'(s - s') (4.6) 
u(s - s') = 5(s - s') (4.7) 

This limit is a simple and useful one, but our main results do not depend on it. 

The propagator J may be evaluated exactly for the case of the simple harmonic oscil- 
lator, V(x) = ^MuP'x 1 '. Introducing X = x + y, £ = x — y, one has 

J(X f ,^ f ,t\X ^ ,0) = F 2 (t) exp(^S-^j (4.8) 



where 



and 



S = K{t)Xfy + K(t)X Q to - L(t)Xo£ f - N(t)X f £o (4.9) 



cj> = A(t)tf + B(t)Z f Zo + C(t)$ (4.10) 

Explicit expressions for the coefficients K, K, L, N, A, B and C are given in Refs. [23,26]. 
F 2 (t) = N/n is a normalization factor, fixed by imposing the condition 

J dxdy 5(x-y) J(x,y,t\xo,y ,0) = 5(xo-yo) (4.11) 

This ensures that Trpt = 1 at all times. On the other hand, tracing over the initial 
arguments of J leads to 

f N 

/ dxodyo S(xq -y ) J(x,y,t\x ,yo,0) = —6{x-y) (4.12) 
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We remark that S is in fact the action of the solution to the boundary value problem for 
the harmonic oscillator with (non-local) dissipation, for which the equation of motion is, 

X + u 2 X + 2[ ds' ri(s-s')X(s') = (4.13) 

Jo 

In the classical limit, we expect that the quantum system reduces to motion described by 
this equation. 

One may also derive an evolution equation for p, for general potentials. Its most general 
form is [26], 

.Bp h 2 fd 2 p d 2 p\ r_ j. . . _ _ / x i 

- ihT(t)(x -y){^-^)- iT(t)h(t){x - yfp 
+ hT(t)f(t)(x-y)(^- + ^) (4.14) 



dx dy / 

Here Vr(x) is the renormalized potential, Vr(x) = V(x) + \M5VL 2 {t)x 2 . The explicit 
forms for the time-dependent coefficients, 5fl(t), T(t), f(t), h(t), are in general rather 
complicated. Explicit expressions for them may be found in Ref. [26] . In the Fokker-Planck 
limit, one has 

9MkT 

r(f) = 7 , M*) = — /(*) = (4.15) 

The first two terms on the right-hand side of (4.14) generate purely unitary evolution 
(but with a renormalized potential) . The third term is the dissipative term, and the fourth 
and fifth terms are diffusive terms. In particular, the fourth term is responsible for the 
process of decoherence discussed elsewhere [5,6,7,8,9]. 



V. TIME EVOLUTION OF I t 



We now study the evolution of I as the density matrix p evolves under the non- 
unitary evolution discussed in the previous section. For simplicity, consider first the unitary 
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evolution of p, without an environment. One has 

K(p,q) = (z\e- iHt p e iHt \z) (5.1) 

where pq is the density matrix at t = 0, and may be pure or mixed. The operators 
e T*Ht^ evo i v j n g pQ forward in time, may be equally regarded as evolving the coherent 
states backwards in time. For a harmonic oscillator, the width a q of the coherent states 
\z) may be chosen to be the width of the ground state (although this choice is by no means 
obligatory). With this choice, the coherent states are preserved under unitary evolution, 
with their centers following the classical evolution: 

e- iHt \p, q ) = \p cl (t),q d (t)) (5.2) 

The same is true for evolution backwards in time, with t — > —t. It is a standard result that 
the transformation from (p, q) to (p c i(t), q c i(t)) is a classical canonical transformation. The 
effect of unitary evolution in (5.1) is therefore to perform a canonical transformation on the 
arguments of fj>t(p, q) at t = 0. It is straightforward to see that our measure of uncertainty 
(2.15) is invariant under canonical transformations of the variables of integration. We 
therefore find that I is constant under unitary evolution for the harmonic oscillator, with 
the above special choice of a q . 

If the width a q is not set to the above special value, then the coherent states are 
not preserved under evolution by the harmonic oscillator Hamiltonian. Likewise for more 
general Hamiltonians. For example, if the initial state is a coherent state, it will spread as 
time evolves, and thus I will increase from its initial value, 1 = 1. Whether / increases or 
decreases, however, depends very much on the initial state. For example, the pure state 
e +iHt \z), which could have a very large value of /, will evolve under e ~ lHt into the coherent 
state | z), possessing the minimum value of /. This "reassembly" of a state sharply peaked 
in phase space from a very spread out state will therefore cause I to decrease with time. 
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One would in fact expect initial states undergoing an initial decrease of I to be just 
as likely as ones undergoing an initial increase, since quantum mechanics is a completely 
time-symmetric theory However, I does in a certain sense capture the intuitive notion 
that "entropy increases", even for pure states, in that it will increase for initial states 
which might reasonably be described as highly organized or special (namely, states that 
are sharply peaked in phase space). 

Now consider the coupling to an environment, as described in the previous section. We 
shall derive an evolution equation for It. We will first use the evolution equation for p, 
(4.13), to derive an evolution for the Wigner function of p, (2.12). Performing the Wigner 
transform of (4.13), one obtains, 

+ E(?)W 2 -^ 

The infinite power series incurred for general potentials makes progress rather difficult. We 
shall therefore restrict attention to the harmonic oscillator, V(q) = \M<J 2 'q 2 , returning at 
the end to a heuristic discussion of the possible effects of more general potentials. Now 
using the expression for p,(p, q), (2.11), one obtains, 

d t JL P d l J , 2u\- d f JL ( a P 2u\ 2 ^u\su\\ 



-& = - £1 + "<4m% -hi- - WWW 



dt Mdq HK ^dp \M uy ' q KJJy J jdpdq 



+ 2T(t)„ + 2T(t) ( p + 4 -)£ + mt)h(t)-g (5.4) 



Here, oo\{t) = u 2 + 5Q 2 (t) is the renormalized frequency. Differentiating the expression 
for 7, (2.15), one obtains, at some length, 

^_ 2r(t) _^_ M4((S2 _ Br(t)/(() )/^I|| 

+ (W ) + 2r ( ^)/^I(|) 2 (5.5) 
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This is the exact result for the time evolution of I for linear systems. 

Now the interesting question is whether we can say anything definite about the mono- 
tonicity properties of /, given Eq.(5.5). First, note that in the case of no environment, and 
for the harmonic oscillator (i.e. ui 7^ 0), it is possible to make the choice 

a l = ^T7~ , °l = \mu% (5.6) 
q 2Muj p 2 y J 

and thus / = 0, as expected. 

The next interesting case to consider is the Fokker-Planck limit, (4.15), in which it is 
again useful to make the choice (5.6), and the second term in (5.5) vanishes. Consider 
the remaining terms in (5.5). The first term is —27 and the coefficient of the last term 
is approximately IM^kT (the a 2 term is negligible in the Fokker-Planck limit). Now the 
question is, what are the relative sizes of the first and last terms in (5.5)? Introduce the 
timescales, 

h 2 1 

tdec= 2M 7 fc7V2' *^ = - (5 ' 7) 
The timescale t^ ec frequently emerges in studies of decoherence, and is therefore called 
the decoherence timescale. We are not of course discussing decoherence per se here, but 
we will use the nomenclature. t re i is the relaxation timescale. On dimensional grounds it 
is clear that the first term will cause I to decrease on a timescale t re i, and the last term 
will cause it to increase on a timescale t^ ec . The important point is that the relaxation 
time is typically very much longer than the decoherence time [6] , so the decoherence term 
will dominate in (5.5). Thus for the harmonic oscillator, with the choice (5.6), and in the 
Fokker-Planck limit, I will increase monotonically for any initial state. 

Now consider the case in which the choice (5.6) is not made. Closely related is the case 
of the free particle, in which u = in (5.5), and a p is arbitrary. The question is whether / 
may be rendered negative by the indefinite term in the integand (which is associated with 
spreading or reassembly). Physically, it is reasonably clear how this may come about. As 
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discussed above, it is possible to choose special initial states that reassemble, at least under 
unitary evolution, and will cause I to decrease. One would expect to be able to identify 
a spreading or reassembly timescale, t s . If the decoherence time scale is much shorter 
than the spreading time scale, one would expect I to increase monotonically, since the 
environment acts before the system has time to undergo reassembly. On the other hand, if 
the spreading time is shorter than the decoherence time, an initial decrease may occur for 
carefully chosen initial states, but this will eventually go over to increase after a time of 
order t^. A similar situation could be expected to hold for more general Hamiltonians.The 
Hamiltonian terms (in (5.3), say) may make / increase or decrease, but eventually the 
diffusive terms will take over and cause / to increase. 

These statements all apply to the high-temperature regime, in which thermal effects 
will eventually dominate. Eq.(5.5) is valid for all regimes, and it would be of interest to 
explore these, although we do not do so here. 

We now have a general picture of the behaviour of / under time evolution. This sets 
the stage for the next Section, in which we derive a lower bound on the behaviour of I. 

Finally, we note that the analogue of Eq.(5.5) for von Neumann entropy is very hard, 
if not impossible, to derive, even for linear systems. Generally it can be obtained only 
if explicit diagonalization of p is possible, e.g., for Gaussian density matrices. For this 
reason, / may be more practically useful than S as a measure of uncertainty, quite simply 
because it is easier to calculate. 



VI. A LOWER BOUND FOR I t 



We now come to the main point of this paper, which is to establish a lower bound over 
all possible initial states for It, thus generalizing (2.16) to include thermal fluctuations in 
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time-evolving non-equilibrium systems. We therefore seek a time-dependent function J™ n 
such that for every time t, 

It > ir n (6-1) 

jrran re p resen t s the least amount of uncertainty the system must suffer, after evolution for 
time t in the presence of an environment. Clearly for consistency, we must have I™ in = 1 
in the absence of an environment. 

To fix ideas, consider first the case of no environment, for which the evolution is 
unitary. The Lieb-Wehrl result is that the information (2.15) at a fixed time is minimized 
by a system in a coherent state, pq = \z')(z'\. A harmonic oscillator initially in a coherent 
state with a width given by Eq.(5.6) evolves so that it remains in a coherent state, and 
therefore It = 1 = I^ flin . It is easy to see that this behavior is very special and cannot be 
realized for other Hamiltonians. This is because Hamiltonian evolution generally does not 
preserve the coherent states. As described in the previous section, for every time r, there 
is an initial state e +iHT \z'), with non-minimal It at t = 0, which evolves to a coherent 
state at time r, there minimizing It. After this, it disperses, and It is no longer minimal. 
It is only minimized at t = r. 

The implication of this is that /™ n is actually an envelope. No particular pq realizes 
the minimum for all time - instead there are a succession of states which achieve the 
minimum. The minimum I™ 1 ™ = 1 is realized, at each time t, by the value of It for the 
initial state e + t \z'); that is, for the initial state obtained by evolving the coherent state 
\z') at time t backwards to t = 0. 

Now consider the situation with an environment, as discussed in the previous section. 
Instead of unitary evolution under e~ lHt , we now have non- unitary evolution under the 
propagator J. As we have seen, interaction with the environment will cause It to increase 
in the long run, but there is the possibility of an initial decrease of It, due to the reassembly 
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effect. We therefore expect I™ m to again be an envelope: there will be many initial states 
which achieve I™ iin for some value of t, but there will be no initial state for which If = I™ m 
for all t. 

To find I™ 171 , we will exploit the Lieb-Wehrl inequality (2.16). It cannot, however, 
be applied immediately to the case at hand. To see why, consider again the case of no 
environment. One is interested in the information (2.15). Application of the inequality 
(2.9) shows that the minimum is achieved for a pure rather than mixed state. One is thus 
minimizing the integral 

7 = _ /tf IW>| 2 1»IW>| 2 (6-2) 
over all square-integrable wavefunctions ifj. The minimum is found to be achieved for 
\ifj) = \z'), a coherent state. If one expresses the state at a later time in terms of unitary 
evolution from its initial value, = e~ iHt \ipo), one has the expression for the information 
at time t 

** = -/§^ |(.|e-^^o)| 2 ln|(.|e-^|Vo)| 2 . (6-3) 

Minimizing this over all square-integrable wavefunctions \i/jq) is easy because e~ t \ipo) 
is itself a square-integrable wavefunction, so the previous result applies, giving = 
e iHt \z'), as discussed above. 

Now we are interested in the more general case in which the propagator is not unitary. 
We would like to know what the new lower bound on the uncertainty is for systems that 
have undergone interaction with the environment for time t. Denoting the coherent state 
density matrix by p z = \z)(z\, and the initial density matrix by po, the information at time 
t is given by 

h = ~ J ^Tr(p z J t (p ))lnTr(p z J t (p )). (6.4) 

Here, Jt(po) denotes the non-unitary evolution of po, Eq.(4.1). For each time t, we seek 
the pq that minimizes (6.4). Differently put, we need to minimize (6.4) over all density 
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matrices of the form pt = Jt(po), where pq is an arbitrary density matrix. The feature 
that distinguishes this case from the Lieb-Wehrl case discussed above is that this class of 
density matrices is smaller than the class of all density matrices, since evolution under 
J is not invertible. It is therefore difficult to characterize the class over which to do the 
minimization. Since Jt(po) is linear in pq, and using the convexity property (2.9), we again 
deduce that the minimizing pq must be pure. This simplifies the problem somewhat, but 
the inconvenience stated still remains. 

To get around this difficulty, we adopt the following strategy. We are interested in the 
quantity, 

ih(p,q) = (Apt\z) 

= J dxdydx dy (z\x)(y\z) J(x, y, t\x , y , 0) po(x ,yo) (6.5) 

where J is the reduced density matrix propagator, pt is then conveniently written in the 
form, 



tH(p,Q) = j dx dy A$(yo,x ) Po(x ,yo) 

= Tr(A z tP0 ) (6.6) 



where 



A t(yo> x o) = J dxd v ( z \ x )(y\ z ) J( x ,y,t\ x oiyoi°) ( 6 - 7 ) 

The quantity A z t is therefore the final density operator \z)(z\ brought back from time t to 
time zero using J. Note, however, that Af is not a physical density matrix, since from 
(4.12), TrA| = ^ (although one does have f ^§A% = 1). Using the Wigner representa- 
tion, (2.12), one may write 

(H(P,q) = J dpdq W A z(p,q) W pa (p,q) (6.8) 

Since J is Gaussian for the linear case considered here, A%(xQ,yQ) and W^(p,q) are also 
Gaussian. 
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Compare this to the Lieb-Wehrl result, (2.16). The latter may be regarded as stating 
that the information of the distribution 

/i(p, q) = 2 J dp'dq' exp ("^p " ^T^j ^ ^ 

is bounded from below by I = 1, with equality if and only if p is a coherent state. Now 
the point is that (6.8) and (6.9) have a very similar form: they are both Wigner functions 
of an arbitrary density matrix with a Gaussian smearing, but the Gaussian factors are 
not the same. Our aim, therefore, is to perform a series of transformations to bring (6.8) 
into the form (6.9), and then apply the Lieb-Wehrl result (2.16). As we shall see, the 
information of \i is not preserved under these transformations, and thus we obtain a non- 
trivial lower bound, different from (2.16), and depending on the quantity A\. The difficulty 
outlined above is avoided because the evolution under J is contained entirely in A%, and 
the minimization is now over all pure po? a well-defined class to which the Lieb-Wehrl 
result may be applied. 

Turn now to the details. Consider first Eq.(6.7). The final density matrix is 

where as in Section IV, X = x + y,£ = x — y. Under evolution backwards in time by the 

non-unitary propagator J it yields, 

i 

A Z t (yo, x ) = j 2 exp (-a (X - 2q ) 2 - + ^0 [r (^ - 2q ) + Pq]J 



where 



(6.11) 



a »=32^A (612) 

9 9 (B-4a 2 NK) 2 

^ = C + 2 ^ JV - 32.IA (6 ' 13) 

^ = ^+4^(4-^1 (614) 
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Here, 



' ' A + ^-\+\k 2 (6.15) 



(These coefficients may be obtained by a straightforward modification of the calculations 
described in Ref.[4]). Also, Pq-,% are the classical evolution of p, q, evolved backwards in 
time under the dissipative equation of motion (4.10). They are given explicitly by 

P0 = ~{NL - KK)q + ^-p (6.16) 

where the various quantities appearing are defined in Section IV. This transformation from 
p, q to Pq, % is non-canonical, because the evolution is dissipative: 

d^Po) = N = ^ 
d(q,p) L 

Performing the Wigner transform, one thus obtains the explicit form of (6.8): 

N fa \^ f /I 2 - 2\ 

xW po (p,q) (6.19) 

We would like to bring this expression into the form (6.9). Introduce 

A = (^) 2 , ix = VSaqifoao) 1 / 4 (6.20) 

Now perform the following canonical transformation on the integration variables, together 
with the same change of variables on pQ, 

q' = M, p' = -ip-2T q) (6.21) 
fi 

Q = M0, P = -(PO - 2F q ) (6.22) 

Eq.(6.19) thus becomes, 
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where we have introduced 

Wp(p', q') = W P0 (w' + 2%', i-) (6.24) 

There arises the question of whether Wp(p', q') defined by (6.24) is still a Wigner function, 
i.e., of whether there exists a density matrix p whose Wigner transform is (6.24). The 
answer is in the affirmative: linear canonical transformations on the arguments of the 
Wigner function are readily shown to correspond to unitary transformations of p. 

The dependence on p, q in the right-hand side of (6.23) resides entirely in p, q, via the 
transformations (6.16), (6.17) and (6.22). It is convenient to write (6.23) as 

N 

mp,Q) = -j-h\PA) (6-25) 

The factor ^ is nothing more than the Jacobean of the transformation from p, q to p, q. 
The transformation (6.21), (6.22) is canonical so the only contribution to the Jacobean 
comes from (6.16), (6.17), whose Jacobean is (6.18). The information of pt, hi is then 
simply related to that of pt, It- It is 

= In (^pj (6.26) 

The distribution pt is almost of the desired form (6.9), but fails to be because of the 
presence of the factor of A. Positivity of the density matrix (6.11), implies that /3q > cxq, 
and in fact equality holds only at t = 0, and thus one has A > 1. One might have thought 
that the next step is to simply scale p' and q' by A2, thus taking A into the Wigner 
function. However, this scaling would lead to a phase space distribution function which 
is not a Wigner function, ie., it is not the Wigner transform of a density matrix. This 
is easy to see: under such a scaling, the degree to which the Wigner function may be 
peaked about a region of phase space becomes enhanced by a factor of A > 1, and thus it 
is possible to violate the uncertainty principle. Wigner functions scaled in this way cannot 
therefore correspond to density matrices. 
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Instead, the next step is carried out using the following simple fact about convolution 
integrals: when two Gaussians with variances o\ and 02 are convoluted, the variance of 
their convolution, 03, satisfies cr| = a\-\-a\. Let us therefore express the Gaussian smearing 
function in fit as the convolution of two Gaussians: 

1 / (p' -P) 2 (q'-q) 2 \ f , , 1 ( (p-p') 2 (q-q'^ 



^X eXP \-^l 2^f) = J dpdq ^ eXP { ^af- ^ 

v 1 ( (p - p) 2 (q-q) 2 \ , R 07 , 

nh(X - 1) 6XP ^ 2(A - 2(A - l)of ) { } 
We may therefore write fit as 



^ ® = ^()bl) / dpdq 6XP [- 2$ -1)4 ~ 2(A - IK 1 ^ Q) (6 - 28) 



where fit is precisely of the form (6.9), with Wigner function Wp(p ,q), given above by 
(6.24). 

The result (6.28) is as close as we can get to casting ji{p, q) in the form (6.9). However, 
the form (6.28) may be exploited: it is the convolution of a Gaussian with the function jit- 
We may therefore appeal to a theorem of Lieb on the information of convolutions [11]. Let 
/ and g be functions defined in L s (R n ), where s > 1, and let f*g denote their convolution. 
Then the information of / * g, I(f * g), satisfies the inequality, 

exp * gfj > exp + exp (^G?)) (6-29) 

Equality holds if / and g are both Gaussians differing only in the location of their centres 
and in an overall scale of their covariance matrices. 

In our case, the Gaussian function in (6.28) has information 

7 = ln(|(A-l)) (6.30) 

Since fit is of the form (6.9), it satisfies the Lieb-Wehrl inequality (2.16), with equality if 
and only if the Wigner function (6.24) is the Wigner function of a coherent state, 

wf /^ 9 r o/A 1 I (p'-p'i) 2 (q'-q[) 2 \ fft ^ 

W po { W + 2-q , -) = - exp — (6.31) 
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Applying (6.29), we therefore have the lower bound on the information of fit, 

h > l + hiQ(A+l)) (6.32) 
Finally, inserting this in (6. 26), we obtain the desired lower bound on If. 



It > 1 + In I 7^7 



i 



a 



(6.33) 



2N 

This is our main result. The right-hand side is the value of I™ 171 at time t. 

Now consider the conditions for equality in (6.33), to determine the initial state which 
meets the envelope at time t. The information of fit achieves its lower bound when (6.31) 
is satisfied, fit is then a Gaussian, differing only from the smearing Gaussian in (6.28) by 
an overall scaling of their covariance matrices. The conditions for equality in (6.29) are 
therefore also satisfied. This means that the inequality (6.33) achieves equality when the 
initial state is given by (6.31). Inverting the Wigner transform, we find that the initial 
state is the pure state, 

1/4 



4(ao/3 ) 1 / 2 
This is a squeezed coherent state. 



*'(•'') = I - U ^ U/ J ex P ( -[2(<*oA))5 + ^r ](x - qi) 2 + % -p lX ) (6.3 I) 



VII. DISCUSSION AND SUMMARY 

We first discuss the properties of the lower bound (6.33). 

Consider Eq.(6.19). We have been seeking the Wigner function W po that minimizes 
the information of (6.19). Loosely speaking, this means finding the Wigner function which 
has the best overlap with the exponential in (6.19), and hence gives the most peaked 
probability distribution nt{p,q)- We found that the initial state doing the job is (6.34), 
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whose Wigner transform is, from (6.31), 

w , s 1 ( (P~ 2r g-M) 2 V? ( ^\ , 71 x 
W*( P , q ) = -e W ^ ^ ^l^-j) ) ™ 

Now consider the exponential function in (6.19). It is the Wigner function of the final 
coherent state evolved backwards by J. The contours of the Wigner function start out as 
circles. Each contour suffers three effects under this non-unitary evolution: it is distorted 
into an ellipse, its axes are rotated, and its area increases. The distortion factor is given 
by [i in (6.20), the amount of rotation is given by Tq, and the area increase is given by A. 
(There is in addition a translation of the contours, but this preserves the information.) 

Now the point is that the Wigner function (7.1) giving the least overlap in (6.19) is 
the Gaussian pure state which matches two out of three of these effects: it has the same 
distortion and rotation factors. It does not have the same expansion factor A - it cannot 
because we know that the minimizing state must be pure, and pure Gaussian states must 
have A = 1. The minimizing state is therefore the state whose Wigner function is close as 
posible to the exponential in (6.19) subject to the constraint that it be pure. 

Turn now to the explicit form of the lower bound. Using the result of Refs. [4,12,26,23], 
it may be shown that in the Fokker-Planck limit, and for short times, one has 

^0 ~i + 2 7 t+ q ' t + 0(t 2 ) (7.2) 
an/ hr 



Setting Oq to the value (5.6), one thus has 



jrmn = 1 + ]n 



(7.3) 



The Fokker-Planck limit involves kT » huR, so 7™ n increases with time. Eq.(7.3) 
indicates that the thermal contributions to the uncertainty principle start to become ap- 
preciable after a time 

- %2 %UR (7.4) 



a%M~/kT 7fcT 
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The important thing to note is that this is the decoherence timescale defined in Eq.(5.7) - 
the timescale on which interference is destroyed by the interaction with the environment. 

Our results should be compared with the work of Hu and Zhang [12]. They calculated 
the usual uncertainty function, (3.10), for the density matrix obtained by evolving an 
initial coherent state for time t in the presence of an environment. They found that for 
short times, and high temperatures, 



It was these authors who first noted, on the basis of this calculation, the significance of 
the decoherence timescale for the comparative sizes of thermal and quantum fluctuations. 
We thus find close agreement with their work. 

This result has a consequence for the decoherence programme. A reasonable question 
to ask in decoherence models is whether there is a regime in which the interaction with 
the environment is sufficient to induce decoherence, yet induces a noise level less than 
that due to intrinsic quantum fluctuations. Our results, and those of Hu and Zhang [12], 
show that this is not the case: in the Fokker-Planck regime, decoherence and thermal 
fluctuations become important on the same timecale. This means, loosely speaking, that 
the uncertainty principle plays little role in these models. 

It is of interest to explore the form of the lower bound in other regimes. Consider for 
example, the low temperature regime. In the Fokker-Planck (high temperature) regime dis- 
cussed above, the diffusion is controlled by the diffusion constant D = 2M^kT. However, 
as argued by Caldeira and Leggett [23], in the low temperature regime the appropriate 
diffusion constant is D = M^Uuju. An order of magnitude estimate on the size of I™ m is 
therefore obtained by substitution of diffusion constants. One thus discovers that in the 
low temperature regime, the environmentally-induced fluctuations (we can no longer call 
them thermal) grow on a timescale 7" 1 , the relaxation timescale. This shows that I™ 171 is 




(7.5) 
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not just a measure of quantum fluctuations of the distinguished system plus thermal fluc- 
tuations of the environment: it also includes the quantum fluctuations of the environment 
(although these are of course negligible in the high temperature regime). 

Another question to ask is whether it is possible to express our new uncertainty prin- 
ciple (6.33) in terms of the usual uncertainty function U. Recall that the Lieb-Wehrl 
inequality (2.16) may be shown to imply the standard uncertainty principle, (1.1), via the 
steps (2.17)-(2.20). Can a similar derivation be carried out in the case of (6.33)? Steps 
analagous to (2.17)-(2.20) can be carried out, and one obtains, 



However, as before this is not the proper form of the uncertainty principle, because the 
variances on the left-hand side also include the variances of the coherent state, Eqs.(2.18), 
(2.19). The final step of minimizing over a q is rather tricky to carry out because the 
right-hand side of (7.6) depends on a q in a non-trivial way, and one ends up with a fifth 
order polynomial in a q . Also, the alternative method suggested by Grabowski [19] cannot 
obviously be generalized so as to apply to this case. Therefore, we do not give an explicit 
form of our uncertainty relation in terms of the variances of p. The possibility of deriving 
such a relation directly (rather than from (7.6)) will be considered elsewhere [28]. 

We should also compare with the work of Paz, Habib and Zurek [7,8], who looked for 
the set of initial states which generated the least amount of physical entropy, S[p], after 
evolution in the presence of an environment for time t. The motivation for doing this is 
that these states are in a sense the ones most stable under evolution in the presence of an 
environment. This is clearly closely related to our work, since we essentially looked for the 
set of initial states with the smallest value of I after time t. Indeed, Paz et al. claimed 
that the minimizing states are coherent states, whereas for us the minimizing states are 
squeezed coherent states. It turns out that the quantity controlling the squeezing, Tq in 




(7.6) 
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(6.34), can go to zero quite quickly (on the timescale u 1 ). In this case we thus see that 
the results are in agreement. 

In summary, we have discussed the properties of an information-theoretic measure of 
uncertainty (1.3) for a class of non-equilibrium quantum systems. Our measure is closely 
related to von Neumann entropy in the thermal regime, but unlike entropy, it supplies a 
non-trivial measure of uncertainty in the quantum regime. It is also easier to work with 
calculationally than entropy. Our main result is the demonstration that, for linear systems, 
our measure has a non-trivial lower bound, the generalization of the uncertainty principle 
to include thermal (or more generally, environmentally-induced) fluctuations for a class of 
non-equilibrium systems. We have examined the form of the lower bound in some regimes 
of interest. A more detailed examination is best carried out numerically, but this is beyond 
the scope of the present work. 
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